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Abstract 

The final state interactions (FSI) model in which soft rescattering 
of low mass intermediate states dominates is suggested. It explains 
why the strong interaction phases are large in the B& — > tttt channel 
and are considerably smaller in the Bj — > pp one. Direct CP asym- 
metries of Bd — > tttt decays which are determined by FSI phases are 
considered as well. 

1 Introduction 

There are three reasons to study FSI in B decays: to predict (or explain) 
the pattern of branching ratios, to study strong interactions, and to forsee in 
what decays direct CPV will be large. In view of this necessity a model for 
FSI in B decays to two light mesons is suggested and explored in the present 
paper. 

The probabilities of three B — > tttt and three B — > pp decays are measured 
now with good accuracy. The C-averaged branching ratios of these decays 
are presented in Table 1 [1]. Let us look at the ratio of the charge averaged 
Bd decay probabilities to the charged and neutral mesons: 

= Br(B d -> p+p-) = Br(B d -> 

p ~ Br(B d - fPfP) ' n ~ Br(B d - tt ti ) ' 1 } 
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Table 1 



Mode 


Br(10- 6 ) 


Mode 


Br(10- 6 ) 


Bd — > n + n~ 

B d -> 7T V 
B u -> 7T+7T 


5.2 ±0.2 

1.3 ±0.2 
5.7 ±0.4 


# d -> P + p- 
B d -> p°p° 
- p+p° 


23.1 ±3.3 
1.16 ±0.46 

18.2 ±3.0 



C-averaged branching ratios of B — > 7T7t and 5 — > pp decays. 



The large difference of i? p and i2„- is due to the difference of FSI phases 
in S — > pp and 5 — > 7T7t decays (see below). In Section 2 we will determine 
the differences of FSI phases of tree amplitudes which describe B — > pp and 
5 — > 7T7T decays into the states with isospins zero and two from the data 
presented in Table 1. As a next step we will suggest a mechanism which 
produces such phases. Once this mechanism is defined it becomes possible to 
calculate FSI phases of decay amplitudes into states with a definite isospin 
(not only their differences). A central question is: what intermediate states 
produce FSI phases in 5-meson decays into two light mesons. In the weak 
decay b — > uu(dd)d in the rest frame of a heavy quark (which is 5-meson 
rest frame as well) three fast light quarks are produced. Their energies are 
of the order of Mb/ 3 and momenta are more or less isotropically oriented. 
The energy of the fourth (spectator) quark is of the order of Aqcd- This four 
quark state transforms mainly into multi pi-meson final state with the average 
pion multiplicity about 9 (this number follows from the experimentally known 
charged particles multiplicity in e + e~ annihilation at E cm = 3GeV multiplied 
by 1.5 * 1.5 in order to take neutral pions and third quark jet into account). 
The total branching ratio of such decays is about 10~ 2 . However such meson 
state does not transform into the state composed from two light mesons 
moving into opposite directions with momenta Mb/2. What meson state 
does transform into two light mesons can be understood from the inverse 
reaction of two light meson scattering at the center of mass energy equal 
to the mass of 5-meson. The produced hadronic state consists of two jets 
of particles moving in opposite directions. Each jet should originate from 
a quark-antiquark pair produced in the weak decay of 6-quark. The square 
of invariant mass of a jet which contains spectator quark does not exceed 
MbAqcd and is much smaller than M B . The energy of this jet is determined 
by that of a companion quark and is about Mb/2. That is why the square 
of invariant mass of the second jet also does not exceed MbAqcd- So for 
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Z?-decays the mass of a hadron cluster which transforms into light meson 
in the final state should not exceed 1.5 GeV. Following these arguments in 
the calculation of the imaginary parts of the decay amplitudes we will take 
into account only two (relatively light) particle intermediate states for which 
branching ratios of 5-meson are maximal. 

In Section 3 we will calculate FSI phases of tree amplitudes describing 
B — > 7T7T decays taking into account pp, irn and nai intermediate states which 
by t(w)-channel exchanges are converted into nn. We will find that large 
probability of B — > p + p~ decay explains about half of FSI phases of B — > nn 
decays. Relatively small probability of B — > 7r + 7r~ decay prevents generation 
of noticeable FSI phase of B — > pp amplitudes through B — > it + it~ — > pp 
chain. 

We will demonstrate that the strong interaction phase of the penguin 
amplitude is opposite to the result of quark loop calculation, which is very 
important for the value of a direct CPV asymmetry C n + n - = C + _ discussed 
in Section 4. Predictions for CPV asymmetries Coo and Sqq will be presented 
in Section 4 as well and the value of the unitarity triangle angle a will be 
extracted from the experimental data on CPV asymmetry S+-. 

Subject of rare B decays is an object of intensive study nowadays and an 
interested reader can find extensive list of references in a recent paper [2] . 



2 Phenomenology; \Sq — 5%\ and |<5q — <^2 



Let us present B — > nn decay amplitudes in the so-called "t- convent ion", in 
which the penguin amplitude with the intermediate c-quark multiplied by 
VubV* d + VcbV* d + V t bV t * d = is subtracted from the decay amplitudes [3] : 



M 



= ^\V ub V: d \mlUf + (0) L-*rJLA 2 e"! + 



+ e'^ ^A e i5 o + 
V6 



v:,v tb 



V ubV* d 



e ipp e i(5" p +5Z) 



Mf 



'd — >7T 7V 



VubV: d 



(2) 



(3) 
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M^-^ = ^|V^V^K/ w / + (0)|^ e -^A 2 e^| , (4) 

where are the elements of CKM matrix, 7 and /3 are the unitarity triangle 
angles and we factor out the product m^f^f+iO) which appears when the 
decay amplitudes are calculated in the factorization approximation. A 2 and 
A are the absolute values of the decay amplitudes into the states with 1 = 2 
and 0, generated by operators 0\ and O2 (tree amplitudes), while P is the 
absolute value of QCD penguin amplitude (generated by operators 3 — 6 
of effective nonleptonic Hamiltonian which describes b quark decays into the 
states without charm and strange quarks). 5%, 5% and 5q are FSI phases of 
these three amplitudes, and it is very important for what follows that all of 
them are different. It is easy to understand why 5$ is different from b\ '■ strong 
interaction depends on the isospin and is different for 1 = and 1 = 2. For 
example, there are definitely quark- ant iquark resonances with 1 = 0, while 
exotic resonances with 1 = 2 should be made from at least four quarks and 
their existence is questionable. The reason why 5$ differs from 5$ is more sub- 
tle. Let us consider the intermediate state made from two charged p-mesons 
which contributes to FSI phases: B d — > p + p~ — > nn. p + p~ intermediate state 
contribution to FSI phases can be large since Br(B d — > p + p~) is big. Both 
tree and penguin induced amplitudes get FSI phases through this chain. Its 
contribution to 5q is proportional to ^{BiB d — > p + p) T / {BiB d — >• 7r+7r~)^ 



t ~ 
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BrB d — > p + p )/(BvBd — > 7r + 7r ) pa 2.1, while that to Sq is proportional to 



(BrB d -> p+p-) P /(BiB d -> n+n-) P . 
How can we determine the penguin contributions to the probabilities 
of B d — > p + p" and B d — > 7r + 7r~-decays? The most straightforward way 
suggested in literature is to extract them from the probabilities of B u — > 
K°*p + and B u — > i^ 7r + decays to which tree amplitudes almost do not 
contribute [4, 5] 1 : 



2 r , 12 

Br(B d -> p+p-) P = 



/p 



A"* 



Ay^ + (l - P f 



T^Br(K°*p + ) w 0.34-10~ 6 

(5) 



1 Contribution of tree amplitudes to these decays comes from the rescattering (B u — > 
if + 7r°)T, -ft'+Tr — » X 7r + , and taking into account CKM suppression of the tree am- 
plitudes of -B — > Kir(K*p) decays relative to the penguin amplitudes we can cautiously 
estimate tree contribution as not more than 10% of penguin one . 
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Br(B d -> 7T+7T-), 




2 

-6 



XJrj 2 + (1 - p) 2 ^Br(ir°7r + ) « 0.59-10 



(6) 

where / p = 209 MeV and f K * = 218 MeV are the vector meson decay 
constants, A = 0.23, rj = 0.34 and p = 0.20 are the CKM matrix parameters 
in Wolfenstein parametrization [6], fx/fi* — 1-2 and the central values of 
Br(5 M -> fT°*p + ) = (9.2±1.5)-10 -6 and Br( J B u -> i^°7r+) = (23.1±1.0)-10 -6 
[1] were used. The accuracy of equations (5) and (6) depends on the accuracy 
of d <-> s interchange symmetry (?7-spin symmetry) of b — > d(s) transition 
amplitudes described by QCD penguin, however when the ratio of (5) to (6) 
is calculated uncertainty factors partially cancel out and we obtain rather 
stable result: instead of being enchanced as in the case of the tree amplitude 
intermediate vector mesons contribution into penguin Bj — > n + 7i~ amplitude 
is suppressed, (Sq) pp ~ 1/2.8(<5J) PP . Taking into account that fraction of 
longitudinally polarized vector mesons produced in B u — > K°*p + decays is 
about 50% we get additional suppression of (<5q) pp by factor \/2. 

Finally, phase Sp comes from the imaginary part of the penguin loop 
with c-quark propagating in it [8]. In order to calculate 5 P let us consider 
corresponding quark diagram. The charm penguin contribution is given by 
the following expression: 



ft'—P^^Mj^d^l/l-^, (7) 

where k is the sum of momenta of two quarks to which gluon radiated from 
penguin decays: k = pi + P2- One of these quarks forms 7r-meson with the 
spectator quark, so neglecting spectator quark momentum in the rest frame 
of 5-meson we have pi = (^n, ^). The second quark forms another 7r-meson 
with d-quark radiated from penguin: p 2 = — ! y-) where < x < 1 is 

the fraction of n + momentum carried by w-quark. Substituting k 2 = xm 2 , 
into (7) and integrating it with the asymptotic quark distribution function 
in 7r-meson ip n (x) = x(l — x) we obtain the value of Sp which depends on the 
ratio iml/m 2 .. In particular, for nib = 5.3 GeV and m c = 1.9 GeV (which 
correspond to the masses of physical states) we obtain Sp « 10°, a small 
positive value. A nonperturbative calculation of Sp described in Section 3 
demonstrates that the sign of Sp can be negative. 

Our next task is to determine the difference of FSI phases Sq — S% (the 
large value of it is responsible for a relatively small value of R n ). If we neglect 
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the penguin contribution, then from (2) - (4) we get the following expression: 



^3 B+ -2B 00 + \^B +Q 
cos(<%-<%) = ^ , 3r+ , (8) 

where B^s are the C-averaged branching ratios, while r /r + = r{Bd) / t(B u ) = 
0.92. Substituting the central values from Table 1 we get \Sq — 5%\ = 48°. 

Penguin contributions to B ik do not interfere with tree ones because a = 
7r — (3 — 7 is almost equal to ir/2. Taking P 2 terms into account with the 
help of (6) (subtracting 0.59 and 0.30 from the first and the second lines of 
Table 1 numbers describing B — > nn data correspondingly) we get: 

\5l-5%\ = 37° ± 10° . (9) 

The accuracy of this 11° decrease of the absolute value of the phases difference 
is determined by the accuracy of (6) and is not high. In recent paper [2] the 
global fit of B — > 7T7T and B — > nK decay data was made. The tree amplitudes 
of B — > 7T7T decays were designated in [2] by T for B — > it + it~ and by C for 
i? — > 7r°7T . According to [2] the difference of FSI phases between C and T 
equals 5 C = -58° ± 10°, \C\ = 0.37 ± 0.05, \T\ = 0.57 ± 0.05 in the units 
of 10 4 eV. The phase shift between the isospin amplitudes is determined by 
these quantities: 

t. fx r \ 3TCsin(-5 c ) 

^-^' zp + rccosfc-c ' (10) 

and substituting the numbers we obtain: 

5 - 5 2 = 40°±7° , (11) 

the result very close to (9). However, the same d <-> s interchange symmetry 
was used in [2] when relating B — > 7T7t and i? — > i^7r decays. Fit [2] was 
made in the same "t-convention" which we use (see the statement at the end 
of page 3 of the paper [2]: "for simplicity, we will assume ... P tc = -Pi«"), 
therefore the obtained results can be directly compared with ours. 

Now let us consider B — > pp decays. According to BABAR and BELLE 
results p mesons produced in B decays are almost entirely longitudinally 
polarized (f L (p +P -) = 0.98 ± 0.03[9], f L (p + Po) = 0.91 ± 0.4 [10], f L (poPo) = 
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0.86 ± 0.12 [11]). For B decays into the longitudinally polarized p-mesons 
we can write formulas analogous to (2) - (4) and we can find FSI phases 
difference with the help of analog of (8). Substituting the central values of 
branching ratios of B — > pp decays from Table 1 we obtain: \Sq — 5% \ = 21°. In 
order to subtract the penguin contribution with the help of (5) we should take 
into account that in B u — > K°*p + decays the fraction of the longitudinally 
polarized vector mesons equals approximately 50% [12], so we should subtract 
0.17 ■ 10~ 6 in case of decay to p + p~ and 0.08 • 10~ 6 for decay into p°p°. In 
this way we obtain: 

\S p - 8 p 2 \ = 20°^ o , (12) 

and the factor 2 difference between (12) and (9) or (11) is responsible for 
the different patterns of B — > pp and B — > im decay probabilities. Let us 
emphasize that while \Sq — 8%\ being only one standard deviation from zero 
can be very small this is not so for | 5q — 5% \ ■ 

3 Calculation of the FSI phases of B — > tttt 
and B — * pp decay amplitudes 

Among three amplitudes of B — > im decays (2)-(4) only two are independent. 
We will calculate FSI phases of B — > 7r + 7r° and B — > it + it~ amplitudes and 
extract from them FSI phases of amplitudes with a definite isospin. 

Our task is to take into account the intermediate state contributions 
into FSI phases. As it was argued in Introduction we should consider only 
two particle intermediate states with positive G-parity to which S-mesons 
have relatively large decay probabilities. Alongside with mi and pp there is 
only one such state: na\. So we will consider pp intermediate state which 
transforms into by n exchange in t-channel, nai intermediate state which 
transforms into 7T7T by p exchange in t-channel and will take into account the 
elastic channel B — > nn — > nn as well. This approach is analogous to the 
FSI consideration performed in paper [13]. However in [13] 2^2 scattering 
amplitudes were considered to be due to elementary particle exchanges in 
t-channel. For vector particles exchanges s-channel partial wave amplitudes 
behave as s J_1 ~ s° and thus do not decrease with energy (decaying meson 
mass). However it is well known that the correct behavior is given by Regge 
theory: g ^ )- 1 . For p-exchange a p (0) » 1/2 and the amplitude decrease 
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with energy as l/\/s. This effect is very spectacular for B — > DD — > nn 
chain with D*(DZ,) exchange in t-channel: a£>*(0) ~ —1 and reggeized D* 
meson exchange is damped as s -2 ~ 10 -3 in comparison with elementary 
D* exchange (see for example [14]). For 7r-exchange, which gives a dominant 
contribution to pp — > Tin transition (see below), in the small t region the pion 
is close to mass shell and its reggeization is not important. 

We will use Feynman diagram approach to calculate FSI phases from 
the triangle diagram with the low mass intermediate states X and Y (see 
Figure 1). Integrating over loop momenta d A k we assume that integrals over 
masses of intermediate states X and Y decrease rapidly with increase of these 
masses. Then choosing z axis in the direction of momenta of the produced 
meson M\ we can transform the integral over ko and k z into the integral over 
the invariant masses of clusters of intermediate particles X and Y 



1 



L B 



dk dk z = / dsxds Y (13) 



and deform integration contours in such a way that only low mass interme- 
diate states contributions are taken into account while the contribution of 
heavy states being small is neglected. In this way we get: 

M^ = M^{8« x 8«Y + iT J x f^) , (14) 

where M^y are the decay matrix elements without FSI interactions and 
TxyLtt-k i s the J = partial wave amplitude of the process XY — > irn (T J = 
(S J — l)/(2i)) which originates from the integral over d 2 k±. 

For real T (14) coincides with the application of the unitarity condition 
for the calculation of the imaginary part of M while for the imaginary T the 
corrections to the real part of M are generated. 

Let us calculate the imaginary parts of B — > nn decay amplitudes which 
originate from B — > pp — > 7T7T chain with the help of unitarity condition 2 : 

/d cos 9 
——M(pp^7Tn)M*(B^pp) , (15) 
327T 

where 9 is the angle between p and n momenta. For small values of 9 or t 
7r-exchange in t-channel dominates and the calculation of Feynman diagram 



2 in this section the phases which originate from CKM matrix elements arc omitted. 
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Figure 1: Diagram which describes FSI in the decay of heavy meson Mq q 
into two light mesons Mi and M 2 . X and Y are the clusters of particles with 
small invariant masses sx,s Y < MqAq C d, k is 4-momentum of a virtual 
particle propagating in t-channel. 



for pp — > 7T7T amplitude with the elementary virtual 7r-meson exchange can 
be trusted, as it was noted above. It was already stressed that p-mesons pro- 
duced in £>-decays are almost entirely longitudinally polarized. That is why 
we will take into account only longitudinal polarization for the intermediate 
p-mesons and amplitudes of -B-decays into tttt and PlPl are simply related 3 : 



/ 18 2 / 23 1 

M B +^ p+p o = - J — , M Bd ^ p+p - = _ Y^2" M Wt • ( 16 ) 

For the amplitude of p + p° — > -k°tt + transition we have: 

<M( P v - -°- + ) = - \ Pl Jc;2_ ml ^p + )^p°) > ( i7 ) 

where p±, k\ and k 2 are p + , ir° and n + momenta. From the width of p-meson 
we get g 2 /lQ'K = 2.85. For the longitudinally polarized p-mesons in their 
center of mass system we have: 

1 



hp + = hp 



2m p 



2 

Tfl 

(t-ml)(l + ^)+m 2 p 



(18) 



where t — (pi — ki) 2 . Changing the integration variable in (15) to t with the 
help of dt = -^-(l — 2jp-)dcos9 and introducing formfactor exp(t/p, ) with 



3 relative negative sign of the amplitudes follows from the expressions for transition 
formfactors in the factorization approximation, see for example [15]. 



9 



the parameter p? ~ 1 GeV 2 we obtain: 



ImM, 



+ 1 



'18.2 



^ 2 2 %2 



J 

/ 



16ttM| * 4m 2 , 



2m 2 

<«-»4><i + ls #>'+ 



2, 2m « 



t — ml 



exp{t/p 2 ) } M B 



(19) 



For p 2 = 2m 2 the contributions of the first two terms in square brackets 
cancel, while the third term gives: 

lmM B ^ n+n o = 
and from (4) we get: 



18.2^ m 2 p 
5.7 16tt4M| 



3.1M B ^ n+7T o , 



(20) 



5UPP) = -4.9° . 



(21) 



Let us note that in the limit Mb — > oo the ratio Br[B d — > pp) / Br(B d — > 7r7r) 
grows as M|, that is why FSI phase 5 2 {pp) (and <5q (pp)) diminishes as 1/M B . 

The analogous consideration of p + p~ intermediate state leads to the pos- 
itive FSI phase of Bj — > 7r + 7r~ amplitude which is enhanced relatively to 
5%(pp) according to (16): 



5l_(pp) = +5.7° , 



(22) 



and for FSI phase of the amplitude with isospin zero in the linear approxi- 
mation we get: 



*5(pp) = K-(pp) + 



A, 



V2A 



(23) 



We are able to extract the ratio A 2 /A from that of C-averaged Bi(B d — > 
7r + 7r _ ), Br(B d — > 7T 7r°) and Br(_B M — > 7r + 7r ), subtracting penguin contribu- 
tion as we did deriving (9): 



at 1 

A 2 



B, + B 



oo 



2 D T 

3^+° 7^ 



1 , 



(24) 
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-f- = 0.80 ±0.09 , (25) 

^2 

and, finally: 

S5(pp) = 15° , 85(pp) - %(pp) = 20° . (26) 

In this way we see that B — > pp — > 7r7r chain generates half of the experi- 
mentally observed FSI phase difference of B — > 7r7r tree amplitudes. 

It is remarkable that FSI phases generated by 5 — > 7T7t — > pp chain are 
damped by Br(£> — > p + p~ , p + p°) /Br(B — > 7r + 7r, 7r + 7r°) ratios and are a few 
degrees: 

^ (7r7r) = ^ * = " L4 ° ' 5 +- (7r7r) = Mi * = L2 ° ' 

(A /A 2 ) pp = l.l, ^(tttt) = 2.9° , ^(tttt) - ^(tttt) » 4° . (27) 

Next we will take into account 7T7t intermediate state. From Regge analy- 
sis of 7T7T elastic scattering we know that good description of the experimental 
data is achieved when the exchanges of pomeron, p and / trajectories in t- 
channel are taken into account [16]. Pomeron exchange dominates in elastic 
7T7T — > 7T7T scattering at high energies. For ctp(0) = 1 the corresponding 
amplitude T is purely imaginary and the phases of matrix elements do not 
change [3]. However taking into account that pomeron is "supercritical", 
op(0) ~ 1.1, we obtain the phase of the amplitude generated by pomeron 
exchange 4 which cancels the phases generated by p and / exchanges for 
1 = 2. For 1 = the sum of p and / exchanges produces the purely imag- 
inary amplitude T and the phase of the amplitude M is due to pomeron 
"supercriticallity" : 

^(tttt) = 5.0° , <S£(7T7r) = 0° . (28) 

In paper [3] the pomeron exchange amplitude was considered as purely 
imaginary. As a result though important for branching ratios phase difference 
5j(7T7r) — 5% {jm) was the same (pomeron contribution being universal cancels 
in the difference of phases) it came mainly from 5J(7T7t) negative value. In 

4 The amplitude of 2 — > 2 process due to supercritical pomeron exchange is T ~ 
(s/s )" p(t) ( 1 + exp{-iirap(t)))/{-sm(Tra P (t))) = (s/s ) (1+A) (i + Att/2), where in the 
last expression t = was substituted and ap(0) = 1 + A was used (A ss 0.1). 
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this way result for the absolute value of direct CP-asymmetry (7^+^- was 
underestimated, see below. 

Finally ira\ intermediate state should be accounted for. Large branching 
ratio of B d -> -decay ( Br(B d -> n ± af) = (40 ± 4) * 10~ 6 ) is partially 
compensated by small pnai coupling constant (it is 1/3 of pmr one). As a 
result the contributions of i\a\ intermediate state (which transforms into tttt 
by p-trajectory exchange in t-channel) to FSI phases equal approximately 
that part of tttt intermediate state contributions which is due to p-trajectory 
exchange. Assuming that the sign of the 7ra\ intermediate state contribution 
into phases is the same as that of elastic channel we obtain: 

^(ttoO = 4° , (^(ttoi) = -2° . (29) 

Summing the imaginary parts of the amplitudes which follow from (21), 
(26), (28) and (29) we finally obtain: 

5^ = 23°, 51 = -7° , = 30° , (30) 

and the accuracy of these numbers is not high, at the level of 50%. 

The analogous consideration of the real parts of the loop corrections to 
B — > tttt decay amplitudes leads to the diminishing of the (real) tree am- 
plitudes by ps 30%, and we can explain the experimentally observed value 
5q — 5% ~ 40° in our model while for pp final state the analogous difference 
is about three times smaller, 8$ — 8% ~ 15°. 

Let us estimate the phase of the penguin amplitude 8p considering charmed 
mesons intermediate states: B — > DD, D*D, DD* , D*D* — > tttt 5 . In Regge 
model all these amplitudes are described at high energies by exchanges of 
.D* (/^"trajectories. An intercept of these exchange-degenerate trajectories 
can be obtained using the method of [17] or from masses of £>*(2007) - 1~ 
and 1^2(2460) - 2 + resonances, assuming linearity of these Regge-trajectories. 
Both methodes give q;d*(0) = —0.8 -. — 1 and the slope a' D * ~ 0.5GeV~ 2 . 

The amplitude of D + D~ — > tt + tt~ reaction in the Regge model proposed 
in papers [18, 19] can be written in the following form: 

T DD ^(s,t) = -| e — (*>r(l - a D *(t))(s/ Scd r D * (t) , (31) 

5 These amplitudes are considered as penguin due to the proper combination of CKM 
matrix elements. 
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where T(x) is the gamma function. 

The t-dependence of Regge-residues is chosen in accord with the dual 
models and is tested for light (u,d,s) quarks [18]. According to [19] s cd 
2.2GeV 2 . 

Note that the sign of the amplitude is fixed by the unitarity in the t- 
channel (close to the D*-resonance). The constant is determined by the 
width of the D* — > Dtt decay: <7q/(167t) = 6.6. Using (14), analog of (15), 
(31) and the branching ratio Br(B -> DD) « 2 • 1(T 4 [20] we obtain the 
imaginary part of P and comparing it with the contribution of P in B — > 
7r + 7r~ decay probability (6) we get S P ~ — 3.5° 6 . A smallness of the phase 
is due to the low intercept of -D*-trajectory. The sign of Sp is negative - 
opposite to the positive sign which was obtained in perturbation theory (7). 

Since DD-decay channel constitutes only 10% of all two-body charm- 
anticharm decays of -B^-meson [20] taking these channels into account we can 
easily get 

5 P 10° , (32) 

which may be very important for the interpretation of the experimental data 
on direct CP asymmetry C + _ discussed in the next section. 



4 CP asymmetries of Bd(Bd) — > tttt decays 

The CP asymmetries are given by : 

r _ 1 - \Ktt\ 2 q _ 2Im(A 7r7r ) _ 2 jf3 ^B^mr /no\ 

7r7r 1 _i_ I \ 12 ' n7T ~ 1 i I \ 12 ' wn ~ M ' ^ > 

where rnr is 7r + 7r~ or tt 7t°. 

From (2) for direct CP asymmetry in Bd(Bd) — ► vr + 7r _ decays we readily 
obtain: 

p 

C + _ = -=sina[V2A sin(5 - 5 - 5 P ) + A 2 sm(5 2 - S - 5p)]/ 

v3 

' ~6~ + 12 + ~Zy/2 C ° ~ ~ V 3 C ° S " C ° S _ ~ ~ 



6 In integration over cos# the region S « 1 dominates. In this region representation 
(31) is valid. 
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where 



A P 

~ cos a cos(<5 2 -6q- 5 P ) + P 2 } , (34) 



VZVtt 



r r - • (35) 

V ubV ud 

In order to make a numerical estimate we should know the ratios A /A 2 
and P/A2. The first one is given by (25) while the second one can be extracted 
from the ratio Br(B u -> K°tt + ) /Bt(B u — > 7T°7r + ) assuming d <-> s invariance 
of the strong interactions: 



Br(5 u ^X°7r+) /iP 2 |^ s T4 6 | 2 



Br(B u -> 7r°7r+) /2|A||K* d K & l 2 
P 



(36) 



0.092(0.009) . (37) 

The numerical values of and are given with good accuracy by factor- 
ization calculation, while P appears to be 2.5 times larger than factorization 
result [3]. In view of this the validity of factor fx in (36) which originates 
from factorization calculation of the penguin amplitude is questinable. If 
factorization of the penguin amplitudes is not assumed then the ratio fx/ fw 
in (36) should be replaced by unity. In this way we get 20% larger value of 
P/A 2 in (37) and we will take this value of uncertainty as an estimate of the 
theoretical accuracy of the determination of P: 

= 0.21(0.04) , (38) 

Taking into account that unitarity triangle angle a ~ 90° and angles 5 and 
5p are of the order of few degrees from (34) we obtain: 

C + _ ~ — 0.28[1.1 sin(5 — 5 — 5 P ) + sin(5 2 — 5 — 5 P )] w 

w -0.56sm((5o + 5 2 )/2 - 5 - 5 P ) . (39) 

In order to determine the lower bound on the value of C + _ let us suppose 
that So = 37°, 82 = 0° (we keep the difference 5 — 5 2 = 37°, as it follows from 
the data on B — > nn decay probabilities (9)), and neglect small values of So 
and 5 P : 

C+_ > -0.18 . (40) 
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Concerning experimental number it could well happen that finally it will 
be considerably below our bound. In this case the result of nonperturbative 
calculation of penguin phase will be confirmed. Substituting in (39) 5 = 
30°, 8 2 = —7° and 5 P from (32) we obtain the following central value: 

C+_ = -0.21 . (41) 

It is instructive to compare the obtained numbers with the value of C H 

which follows from the asymmetry Acp{K + tt~) if d <-> s symmetry is sup- 
posed [21]: 



C+- = MH A CP (K+n 



fx) ' T(B—>n + n ) sin(7) 



V t 



td 



= 1.2(~ 2 )(-0.093 ± 0.015) 19 ' 8 Sm82 ° 0.87 = -0.24 ± 0.04 . (42) 
v ; 5.2 sin 60° v ; 

Let us note that one factor f n / fx in the last equation appears from the 
matrix element of the tree operator, the second one - from the matrix element 
of the penguin operator. If because of nonfactorization of penguin amplitudes 
we will omit the factor which appears from the penguin [5] , then the numbers 
in the right-hand sides of (40, 41) and (42) will become 20% smaller. 

The experimental results obtained by Belle [22] and BABAR [23] are 
contradictory 

C Beiie = _o.55( . 09 ) c babar^ _ .21(0.09), (43) 

Belle number being far below (40) and (41). 

For direct CP asymmetry in B d (B d ) — > 7r°7r° decay from (3) we readily 
obtain: 

C 00 = -^P sin a[A sm(5 - 5 - 5 P ) - V2A 2 sm(5 2 - S - S P )]/ 
A 2 A 2 \Fl Pi 

I [ + -j- - y a ° A2 cos ^° - ^) - y 2 A ° P cos a cos ^° ~ 6 °~ 6p "> + 

2 ~ 
+ -^A 2 P cos a cos(5 2 - 5 - 5 P ) + P 2 } , (44) 
V3 

Coo ~ -1.06[0.8 sin(5 - S - 5 P ) - 1.4 sin(<5 2 - 6 - 6 P )\ w -0.6 , (45) 
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considerably smaller than C + _. This unusually large direct CPV (measured 
by | Coo |) is intriguing task for future measurements since the present exper- 
imental error is too big: 



C™ per = -0.36(0.32) . (46) 

Belle and BABAR agree now on the value of another CPV asymmetry 
measured in B d (B d ) -> tt+tt- decays: S+E er = -0.62 ± 0.09 [22, 23]. From 
this measurement the value of unitarity triangle angle a can be extracted. 
Neglecting the penguin contribution we get: 

sin2a T = S + _ , (47) 

a T = 109° ± 3° . (48) 



Penguin shifts the value of a. The accurate formula looks like: 

Al A\ A A 2 
6 12 3^2 



[sin2a(^ + # + 4 2 #^s(5 -5 2 ))- 



A 2 P ~ 2 

— -j= sin a cos(<52 — So — 5p) — y g^o-P sin a cos(5o — — &p)\I 

' ~6~ + 12 + 1^/2 C ° S _ _ V 3 C ° S a ~ ~ ~ 

A P 

- —=r cos a cos(5 2 - 5 - 5p) + P 2 ] , (49) 

v3 

and since all the phase shifts are not big the values of cosines in (49) are 
rather stable relative to their variations. For numerical estimates we take 
5 = 30°, 5 2 = —7° and neglect 5 and Sp. In this way we get: 

(a)^ = 88° ± 4°(exper) ± 5°(theor) , (50) 

where the first error comes from uncertainty in S e + per while the second one 
comes from that in the value of penguin amplitude, (38). Relatively large 
theoretical uncertainty in the value of P does not prevent to determine a 
with good precision. 

The relative smallness of penguin contribution to B — > pp decay am- 
plitudes allow us to determine a with better theoretical accuracy from the 



16 



experimental measurement of (S^ ) pp just as it was done in [24]. With the 

help of (5) we obtain: 

(^-) pp = 0.060(0.012) , (51) 

A 2 

where the same 20% uncertainty in extracting penguin amplitude is supposed. 
Using the ratio (A /A 2 ) pp determined in (27) from the (49) neglecting strong 
phases (which are much smaller than in the case of B — > mi decays) and 
taking into account the recent experimental result (S^ er ) pp = —0.06 ± 0.18 
[1] we obtain: 

{a) pp = 87° ± 5°(exper) ± l°(theor) . (52) 

Let us point out that considerably larger theoretical error quoted in [4] 
follows from the larger theoretical uncertainty in the value of penguin ampli- 
tude assumed in that paper. 

Our results for a should be compared with the numbers which follow from 
the global fit of unitarity triangle [6, 7]: 



CKM fitter = (QQ^^) , = (93 ± 4)° . (53) 



a 

We conclude this section with the prediction for the value of CPV asym- 
metry S q: 

Al A 2 2 V2A A 2 
S 00 = [sin 2a(-^ + — cos(5 - b 2 )) + 

2A 2 P ~ l~2 

+ -j= sin a cos(5 2 — bo — bp) — y s m ol cos(5 — b — bp)}/ 

. ,A 2 A\ V2A A 2 .. . , [2 . ~ ~ . , 

/ [— + — cos(d - b 2 ) - W-A)Pcosacos(d - o - 6 P ) + 

2AoP 

+ -^cosacos(5 2 -So -bp) + P 2 ] = 0.70±0.15 , (54) 
v3 

a large asymmetry with the sign opposite to that of S+_. 
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5 Conclusions 



FSI appeared to be very important in B — > tttt decays. 

The description of these interactions presented in the paper allows to 
explain the experimentally observed difference of the ratios of decay proba- 
bilities to the neutral and charged modes in B — > nn and B — > pp decays. 

Rather large absolute value of direct CP asymmetry C + _ (if confirmed 
experimentally) will be a manifestation of the negative sign of penguin FSI 
phase in accord with nonperturbative calculation and opposite to perturba- 
tive result. 
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